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We develop a microscopic and gauge-invariant theory for collective modes resulting from the phase 
of the superconducting order parameter in non-centrosymmetric superconductors. Considering var¬ 
ious crystal symmetries we derive the corresponding gauge mode caG(q) and find, in particular, 
new Leggett modes uJL(q) with characteristic properties that are unique to non-centrosymmetric 
superconductors. We calculate their mass and dispersion that reflect the underlying spin-orbit cou¬ 
pling and thus the balance between triplet and singlet superconductivity occurring simultaneously. 
Finally, we demonstrate the role of the Anderson-Higgs mechanism: while the long-range Coulomb 
interaction shifts WG(q) to the condensate plasma mode a;p(q), it leaves the mass Ao of the new 
Leggett mode unaffected and only slightly modifies its dispersion. 

PACS numbers: 74.20.-z, 74.70.-b, 71.45.-d, 74.25.N- 


Introduction. Owing to the Pauli exclusion prin¬ 
ciple in single-band superconductors spin-singiet 
(even parity) and triplet (odd parity) pairing cor¬ 
relations never occur simultaneously. Important ex¬ 
amples are spin-triplet odd-parity pairing correla¬ 
tions in superfluid ^He [i|, U, triplet superconduc¬ 
tivity in Sr2Ru04 Q, as well as unconventional 
singlet pairing correlations in heavy Fermion sys¬ 
tems [5 and cuprates @ . A necessary prerequisite for 
a clear singlet-triplet distinction is, however, the ex¬ 
istence of an inversion center. The discovery of the 
bulk superconductors GePtaSi (tetragonal [g) and 
Li2Pd2;Pt3_a;B (cubic Q), without inversion symme¬ 
try, to give only two examples, has therefore initiated 
extensive theoretical and experimental studies. The 
Rashba-type spin-orbit coupling caused by the ab¬ 
sence of an inversion center implies (i) the lifting of 
the band degeneracy associated with a splitting into a 
two-band structure and (ii) the superposition of both 
singlet and triplet contributions to the superconduct¬ 
ing gap [i,i]. 

The breaking of a continuous symmetry in super¬ 
conductors is associated with the occurrence of a 
gauge mode which is necessary to restore the charge 
conservation. Furthermore, in analogy to the Joseph- 
son effect, Leggett predicted the appearance of a new 
collective excitation in s-wave two-band supercon¬ 
ductors, which corresponds to an out-of-phase os¬ 
cillation mode of the phase difference of the cou¬ 
pled condensates [l^. So far, the Leggett mode has 
been only observed in MgB2 ll|, but several pre¬ 


dictions for other s-wave superconductors have been 
made [l2l - [l3| . In non-centrosymmetric superconduc¬ 
tors (NGS), however, where a complex mixing of sin¬ 


glet and triplet superconductivity occurs, it is not a 
priori clear whether a Leggett mode exists Q. 

In this letter we use a microscopic theory to demon¬ 
strate the existence of Leggett modes in NGS. For 
this purpose, we calculate all order parameter collec¬ 
tive modes associated with the condensate phase dy¬ 
namics. For the first time we provide analytic expres¬ 
sions and numerical calculations for the gauge mode 
WG, the mass and the dispersion of Leggett’s collec¬ 
tive mode a;L(q), as well as for the frequency wp of 
the condensate plasma mode. The interesting inter¬ 
play of these collective modes is studied in connection 
with the electromagnetic response of the pair conden¬ 
sate, with special emphasis on the participation of the 
collective modes in the Anderson-Higgs mechanism 
iSH- We emphasize the calculation of the mass Aq 
of various Leggett modes that depend strongly on the 
singlet-to-triplet ratio and may be observable by Ra¬ 
man or Brillouin scattering experiments. 

Model description of NCS. The Hamiltonian for 
noninteracting electrons in a non-centrosymmetric 
crystal reads 


'fl — ^ ^ ^k(7 “t” Tk ' 
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where ^k represents the bare band dispersion, cr, a' = 
t, i label the spin state and t are the Pauli ma¬ 
trices. The second term describes an antisymmetric 
spin-orbit coupling (ASOG) through the vector 7^. 
In NGS two important classes of ASOGs are realized 
which reflect the underlying point group Q of the crys¬ 
tal. We shall particularly be interested in the tetrag¬ 
onal point group (74^ (relevant for GePtaSi) and the 
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cubic point group 0(432) (applicable to the system 
Li 2 PdxPt 3 _xB). For Q = the ASOC reads 

7k = 7-L(k X e^) + . (2) 

In the purely 2D case (yn = 0) one recovers the Rashba 
interaction. For the cubic point group Q = 0(432) 
reads y^ = yik - + kDe^ + ky{kl + kl)ey + 

kz{kl + kl)e^] . 

What are the consequences of the ASOC? First, 
diagonalizing the Hamiltonian, one finds the energy 
eigenvalues = Ck + Mll7kll with /i = ±1 which 
correspond to a lifting of the band degeneracy be¬ 
tween the two spin states at a given momentum hk.. 
This band splitting is responsible for the two-band 
structure characteristic of NCS metals. Second, the 
presence of an ASOC invalidates the classification of 
the superconducting order parameter with respect to 
spin singlet (even parity) and spin triplet (odd par¬ 
ity). Thus, in general, a linear combination of the gap 
on both bands is possible. Sigrist and co-workers have 
shown that most likely y^. orientates parallel to the 
d-vector of the triplet part Q. Thus, we can simply 
write the gap function on the two bands in terms of a 
singlet (As) and a triplet (A^^) amplitude: 


Ak^ = As(T) + /iAi,(r)/k , (3) 

with /k = ||ykll/[(ll7k'lP)Fs]^^^ - 0. ’^ere (.. .)pg 
denotes the Fermi surface (FS) average [l7|. Thus, in 
short, while for all superconductors having an inver¬ 
sion center either singlet or triplet pairing is realized, 
in NCS singlet and triplet pairing occurs simultane¬ 
ously. Simply speaking, the resulting ASOC may drive 
e.g. s- plus p-wave pairing on one band while s- mi¬ 
nus p-wave is established on the other, leading to new 
collective modes. 

Nonequilibrium Kinetic Theory for NCS. In order 
to calculate the dynamical properties of NCS we con¬ 
sider the response to a scalar electromagnetic poten¬ 
tial In addition there contributes a charge 

fluctuation term, which accounts for the action of the 
3D long-range Coulomb interaction Vq = 47re^/q^ 
within the RPA, i.e. x = where 

X is a generalized response function. Then, the re¬ 
sponse to the perturbation Sf = e(/)(q, w) -I- V^5n(q, w) 
with Sn being the total density response of the sys¬ 
tem, is described by a generalized momentum dis¬ 
tribution function Upp, which is a 2 x 2-matrix in 
Nambu-space using again the band basis of Eq. m 
with p = ±1. At the same time the perturbation 6(^ in¬ 
duces fluctuations 



Figure 1. (color online) Illustration of various calcu¬ 
lated collective modes (T=0) common to all NCS. The 
Anderson-Higgs mechanism shifts the gauge mode ljg 
(dashed line) to the plasma mode cjp usually lying in the 
pair-breaking continuum. The new Leggett modes (solid 
green lines) unique to NCS are only slightly changed by 
this process (not visible) and the mass Aq remains un¬ 
changed. Importantly, in some cases Aq —>■ 0 is possible, 
see discussion of Fig. [3] thus, the Leggett modes might be 
easy observable. Note that the slope of the Leggett modes 
depend on the ratio t = Atr/Ns as discussed in connection 
with Fig. [2] 


of the pairing amplitude pk/x, as well as the important 
phase fluctuations of the superconducting order pa¬ 
rameter = i[^Ak^jJ^ - which 

we will later use to determine all collective modes. The 
Fourier transformation of describes the evolution 
of the system in space and time after perturbation df. 
However, it is convenient to stay in (q, a;)-space and 
solve the von Neumann equation [l^ 






1=0 

pp"’ip"p'l 


( 4 ) 


in the clean limit, where p = 7i(k4-q/2), p' = 
?i (k — q/2) and the 2x2 energy matrix ^p„p, have 
been introduced. The simplest way to solve Eq. o is 
to make the following ansatz: 


ZIppx = 21k/x(q>w) = ngySq^o + 6nky{q,uj) 

= ilAo + ^i,^(q,a.) 


with the nonequilibrium quantities 




Suky ySgky \ .. 

-Sn-uy ) ’ -kfX 


5C, ySAuy \ 
-Sf ) 


After some lengthy, but straightforward calculations 
[supplement material, Eqs. (A.7)-(A.ll)] we obtain 
from the off-diagonal components of Eq. (|3]) the rela- 
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tion between fluctuations of the pairing amplitude and 
fluctuations of the superconducting order parameter: 

2 2 ( 6 ) 
- (q ■ VkA.) lAk^—— . 

Here, we have identified the condensate response 
function 

X _.a 2 - (q ■ Vk^)^] + ’^‘kM(q ■ VkA.)^ , S 

‘''^(q-Vk,) = [c.2-4^£J-a22[c.2-4£2J (0 

with Vk^ = d^i^f./dh'k, = tanh(£;k^/2fcBr)/2£'k^, 
Ekfi = and ^k,. = with 

momentum distribution function rik^. An important 
property of the condensate response is the sum rule, 
which generates the condensate density -^p/i = 
^0 Y.^L (^pm)fs = ^oA , with No = NyI2 being the 
DoS for one spin projection. As we will show in supple¬ 
ment material [Eqs. (A.16)-(A.18)] the total particle 
density 5n obeys the conservation law ui5n — q • j = 0 
only, if all phase fluctuation modes of the order pa¬ 
rameter in Eq. m are properly accounted for. 

Finally, we find from the diagonal components of 


Eq. (H]) the density response of NCS: 




/ (q-Vk^jVkM 

_ (q . 





2Ak^ 


( 8 ) 


with ipkfi. = ‘f’k/i — Ak/i being the quasiparticle re¬ 
sponse. Since we are only interested in the response 
of the superconducting condensate Sus, we may ig¬ 
nore quasiparticle contributions oc ipkfj, in Eq. ([5]). 
Then, the density response function simplifies to 
Srikfx = -Ak^i^C + a;Ak^i^A[.“V2Ak;j. Hence, the 
condensate density response Sug = X^k/x '^n-k/j is ex¬ 
clusively determined by Ak/x- In other words, we find 
that the frequency- and wave-vector dependence of 
(5ns(q, w) contains all information on the relevant or¬ 
der parameter collective modes in NCS. Finally, com¬ 
bining Eqs. (IH) with both the superconducting gap 
equation Ak^ = variation 


( 9 ) 

piy 

(with rj(p being the pairing interaction [l^) leads to 
the main result of our analysis (q = q/|q|): 


(5ns(q, w) = A^oA- 


a;^(q)[a;^-<(q)] 


i(q) -h w^(q) -h wg(q)] -h [a;|,(q) -h w|(q)] uj^{q) 


3(/i(q,w) 


( 10 ) 


New collective modes. From the denominator of 
Eq. (na we can draw important conclusions which 
are summarized in Fig. [T] In analogy to neutral sys¬ 
tems we first consider wp(q) —>■ 0 and find two poles 

wi = a;G(q) + O gauge mode 

W 2 = w£(q)-I-C Leggett mode 

with WG(q) being the characteristic gauge mode of 
NCS with w^(q) = E^(Ap/x(q-Vp/x)2)pg/A . Fur¬ 
thermore, we discover the Anderson-Higgs mechanism 
for the gauge mode in NCS shifting it to the plasma 
frequency, i.e. a;p(q) = u;|,(q) -I- WQ(q). Thus, after 
Coulomb renormalization, we find: 

u)l = a;p(q) -I- a;G(q) -f C* f ‘*' 2 ) plasma mode 

; 2 , (12) 

ujI = Wl (q) -I- O Leggett mode 

with wp(q) being the characteristic conden¬ 
sate plasma frequency of NCS with a;p(q) = 
47rne^ 3 (Ap^(q • p)^)pg/m. It is important to 
note that the full condensate density response Sug as 


described by Eq. is also manifested in the con¬ 
densate dielectric function e = 1 — VqSni^'^ /e(f>^ with 
= Sns{ujp{q) —> 0). All in all, our new results for 
the gauge mode and plasma frequency generalizes the 
known solutions for ordinary two-band superconduc¬ 
tors which can be obtained in the limit /k = 1 [^ . 

The second pole in Eq. (ITOl) leads with u;p(q) —>■ 0 
to Eq. (Hill determining the new Leggett’s collective 
modes WL(q) in NCS corresponding to oscillations in 
the relative phase of the superconducting condensates. 
The exact analytical result for ^^(q) is too lengthy 
to be shown here and thus can be found in the sup¬ 
plement material [see Eqs. (B3)-(B5)]. Instead, we il¬ 
lustrate its dispersion (for different t = At^/As) in 
Fig. [Hand calculate its slope (as an example for C' 4 „) 
in Fig.jH As expected, we find for all point groups con¬ 
sidered the dispersion (^^(q) — Aq) oc jqp. The slope, 
however, depends on the ratio t = Atr/^s- Thus, in 
Fig. UK a) we show the slope of the Leggett mode exem- 
plarily for the tetragonal point group [see Eq. ([2])] 
along the q^;- and qj,-direction. The calculated up- 
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Figure 2. (color online) Slope of the dispersion of the 
Leggett mode for NCS systems with Civ point group sym¬ 
metry as a function of the unit vectors (ix,ciy,<lz- (a) 
Comparison of the slope for t = 0.5 (upward parabola) 
and t = 1.5 (downward parabola), (b) slope along the qz- 
direction for various t — Atr/Os¬ 


wald parabola corresponds to t = 0.5 while the down¬ 
ward parabola corresponds to t = 1.5, respectively. 
For t = 1 one finds a constant slope of 1/3 (indepen¬ 
dent of Qa; and Qy, not shown). The three resulting 
dispersions are illustrated schematically in Fig. [TJ In 
Fig.EJb) we show the slope along the qz-direction for 
various t which reveals a non-monotonic behavior for 
fixed q^. In contrast, for the cubic point group 0(432) 
we find in all directions ^^(q) — Aq = |vp|qp inde¬ 
pendent of t (not shown), since the underlying ASOC 
is isotropic to leading order. This would correspond to 
the curve with t = 1 in Fig. [T] 

From (q) in Eq. dill) , we find the mass Aq of the 
Leggett mode 

Aq = = 0) = 47ncs^sAtrT 7 ^ , (13) 

A0A2 — Ai 

where the definitions 

A„ = ^ (Ap^(A./p)”/A^^)ps (14) 

/i —±1 



Figure 3. (color online) Normalized mass of the Leggett 
mode for NCS systems with Civ point group symmetry 
for fixed As = 0.1 as a function of the mixing term and 
various \tr'- \tT = 0.025 (upper solid line), Atr = 0.05 
(middle solid line) and Atr = 0.075 (lower solid line). The 
dashed lines correspond to Eq. m which is an analytical 
solution in the limit of small t. 


have been used. Here, 7„cs represents the coupling 
strength of the Leggett mode, which we will calcu¬ 
late below. In order to determine Aq we need the 
exact solution of the coupled self-consistency equa¬ 
tions of the superconducting gap functions [see also 
Eq. ([3])]: Ak^ = X/pi/=±i with gp,^ = —dp^^Ap,^ 

being the pairing amplitude and Opy has been de¬ 
fined together with Eq. 0. We choose the generalized 
two-gap weak-coupling pairing interaction of Ref. 21 j 


rC; = -{r. 


r rn (M/k + 


0(eo - |?k/i|)0(eo - l^pi^l) and obtain 


-A-i 


A, 

Atr 


/p)} 


with Aq, = N{0)Ta,a = s,tr,m, 

; A-i = 


A = 


A« 

Arr 


Am 

Air 


1 


Atr Ar 
Am, A e 


(15) 


(16) 


and (a^/p)")fs- Aote that one obtains 

the ordinary two-band case if Si —)■ 0 22|. Equa¬ 
tions dlSI-dlll) have the advantage that the exact re¬ 
lation 


7ncs — Am/|A| + ill (17) 

holds and thus determines the coupling constant in 
Eq. (fTSl) . Thus, for given A^, Atr, Am a numerical ex¬ 
act solution of Eq. (fT51) is always possible: the result¬ 
ing exact gap function Ak^ needs to be inserted in 
Eqs. (US and (fT31) to determine Aq [2^ . 

In Fig. |3] we show results for the Leggett mass Aq 
for fixed A^ = 0.1 as a function of Am- While for a 
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small triplet contribution (upper solid line) Aq in¬ 
creases monotonically, we find a non-monotonic be¬ 
havior of the mass for increasing (middle solid 
line). Finally, if ~ Xtr we obtain the important 
case that Aq can become zero (lower solid line). Phys¬ 
ically, this corresponds to a partly vanishing gap on 
one of the Fermi surfaces [see Eq. © and Ref. |24|]. 
Also displayed in Fig. |3] is the analytical solution in 
the limit of small t (dashed lines) 

Ag = 2A; [1_ (3 ,18, 

This might help experimentalists to estimate in which 
materials the new Leggett modes are most easiest ob¬ 
servable. 

Finally, we return to the Anderson-Higgs mecha¬ 
nism. What is its role for the new Leggett modes? 
First, we conclude that the Leggett mass Aq is un¬ 
changed, since the r.h.s. of Eq. (IT^ does not depend 
on WG- Physically, this corresponds to the fact that 
the Meissner effect in the presence of a new Leggett 
mode is unchanged. Second, we find that the disper¬ 
sion of the Leggett mode is only slightly changed. To 
see this, one needs to consider the difference in wf be¬ 
tween Eqs. m and m- Since a;G(q —0) —)■ 0 and 
wp 3> wl, the higher order corrections nearly vanish. 
The resulting is also very small [see supple¬ 

ment material Eq. (C.4)]. Thus, we conclude that the 
dispersion of the Leggett mode and the results shown 
in Fig. [2] are nearly unchanged due to the Anderson- 
Higgs mechanism [25j. 

In conclusion, using a gauge-invariant theory of 
superconducting phase fluctuations in NCS we have 
demonstrated the existence of Leggett modes and cal¬ 
culated their characteristic mass and dispersion for 
various crystal symmetries. Both properties reflect the 
underlying spin-orbit coupling and depend strongly 
on the singlet-to-triplet ratio. Furthermore, we have 
calculated the corresponding gauge modes and clar¬ 
ified the role of the Anderson-Higgs mechanism for 
collective modes in NCS. 
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A. KINETIC THEORY FOR NCS. 


Model description in equilibrium. A non-centrosymmetric superconductor (NCS) is described in equi¬ 
librium by the Hamiltonian which is given by Eq. (1) in the main text. In order to include the pairing 
correlations into the description, we extend Eq. (1) to include the gap matrix as an off-diagonal element 
of an energy matrix in Nambu space. In the presence of an antisymmetric spin-orbit coupling (ASOC), 
represented by the vector the 4x4 energy matrix has the following form in the spin representation: 


r = 

Ek 


Ckl + 7k • T 


-K-kl + 7-k ■'t] 


(A.1) 


In order to account for the two-band structure occurring in NCS systems in the limit of large spin-orbit coupling, 
it is convenient to perform a unitary transformation of into the helicity-band basis or simply band basis. The 
transformation from spin to band basis is described by the matrix Uk, which has the property 

Ul(7k-r)Uk = ||7kl|r^ (A.2) 


and which is obtained in the form of a SU(2) rotation 


Uk 



■T 


COS 0-y = 7^ • z ; n.y = 


7k X z 

Il7k X z|| 


(A.3) 


that corresponds to a rotation in spin space into the z-direction about the polar angle 0.y between 7 ]^ and 
z. Here, r denotes the vector of Pauli spin matrices. A straightforward extension of this transformation into 
Nambu space reads [S.l] 




/ a+ 

0 

0 

Ak+ 

\ 

0 

a- 

—Ak- 

0 

0 


-a- 

0 



0 

0 

-a+ 

/ 


Uk = 


Uk 0 \ 

0 u*; 


(A.4) 


with the energy values ^k^i = Ck + /.t||7kll and the gap functions Ak^ = As{T) + pAtr{T)fi^ also given by Eq. (3) 
in the main text. Introducing a band-index /i = ±1, one may write the equilibrium energy matrix in the band 
basis in the compact form: 


tO _ f ?k/x /^Ak/i \ 

-k/. “ pAl^ -^_k^ ) 

In analogy, one can find for the equilibrium density matrix: 

0 f 2 M^k/i^k/i A 


(A.5) 


(A.6) 


Nonequilibrium Kinetic Equations. The action of an external perturbation S( = e^(q, to) + Pq(5n(q, w) 
leads to the deviation of the density matrix, as well as the energy matrix, from its equilibrium value. An NCS is 
now described in the band basis by a generalized momentum distribution function u)) , and an energy matrix 
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,, respectively. A collisionless quantum dynamics is given by the von Neumann equation [see Eq. (4) in the 
main text]: 




El. 


,] = o 

’ 2.p"p'j 


(A.7) 


This equation can be linearized by using the ansatz of Eq. (5) given in the main text. This leads to 




•k/j,” 


Ik/.A 


(A.8) 


with the equilibrium quasiparticle energy and the distribution function defined in Eqs. (IA.5I) and (IA.6I) . 
respectively. The momentum and frequency-dependent deviation from equilibrium can be defined in the appro¬ 
priate way as 2 X 2 matrices in the Nambu space: 


Sllkt, = 


f ^^k/i 

V -<5n_k^ 


and = f 


(A.9) 


with 5^k/i = = <^C- Thus, the equation (IA.8I) represents a set of eight equations in the band basis [S.2] 

(with the band index /i = ±1). Furthermore, it is convenient to decompose the diagonal elements of the energy 
and density deviation matrices according to their parity with respect to k ^ — k 


^ + S(5n_k^) 

f(s) 


(A.IO) 




with the labeling s = ±1. By analogy, the off-diagonal components are decomposed into their real and imaginary 
parts: 


^9kl = ^ ( <^5k, 


Jk^i 2 

= - 
k/i 2 


A. 


k/i 


SAuij 


l^k/x 

A 








(A.ll) 


where represents the amplitude fluctuations and the phase fluctuations of the order parameter. 

After these specifications the off-diagonal components of the Eq. (IA.8|) simplify to [S.2]: 


= -1 


^9kJ + ^ S( [w (q • Vk^) ]Ak/i ^^2 


A- 


k/i 


- (Vk/. • q)^ - 4A2 


k/i 




^k/i 


wAk, 




ik/i 


k/i 


whereas for the diagonal elements one gets: 


= 


= 


2A: 
(5A|“^ 

— - + (vk^ • q) Ak/^-v-r 

w^-(q-Vk^) ^^k/i 


- (q- Vk^x) 

(-) _ w (q • Vk^) (fikfj. 


SA 


(+) 

k/i 


(q ■ Vk^)^ </Jk/x , Ea , ^ 

- - -% - Ak^ oC + t^Ak^- — 


k/i 


(A.12) 

(A.13) 

(A.14) 

(A.15) 


Equation (|A.13p describes the important relation between fluctuation of the pairing amplitude and the phase 
fluctuations of the superconducting order parameter [see Eq. (6) in the main text]. The density response to a 























(+)l 


scalar perturbation (5C is given by Eq. (I A. 141) [corresponding to Eq. (8) in the main text with drik/^ = 


Conservation law. One strength of the matrix kinetic equation approach lies in the straightforward physical 
interpretation of its results. In addition, the gauge invariance of the whole theory can be demonstrated easily 
if all phase fluctuation modes of the order parameter are properly taken into account: As one can see from 
Eqs. (|A.14I) - (IA.15I) the density distribution functions are directly connected with the phase fluctuations 

of the order parameter The combination of the results from Eqs. (IA.14|) - (IA.15I) yields together with the 

subsequent integration over the momentum space k to the continuity equation 


wdn - q-j = ^ Ap^ < w^-(q-Vp^) 


6A 


(-) 

PH 


J 2Ar 


-wdC 


(A.16) 


which at first glance displays a non-vanishing right-hand side. However, by using Eqs. (IA.12I) - (IA.13|) and the 
variation of the energy gap equation 




-) 

pi/ 


(A.17) 


pv 


one finds after a straightforward, but lengthy calculation: 


uj5n — q • j = 0 . 


(A.18) 


Thus, the particle conservation and, associated with it, the gauge invariance of the theory are satisfied within 
the framework of the matrix kinetic theory. 

B. NEW COLLECTIVE MODES 

The collective excitations of a non-centrosymmetric system can be obtained from the condition, that the 
denominator of Eq. (10) vanishes, i.e. 


- [a;^(q) -f w^(q) -f w2(q)] -f [u;|(q) -f a;^(q)] WL^(q) = 0 


(B.l) 


with WG(q) = S/i • Vp/x)^)pg /A and w|,(q) = 47rne^ 3 ' p)^)fs 


abbreviation 


and 


UJ- 


, 1 , \ _ a 2 , Q;oQ2a^qOa^q2 - gj^ql 

^ ^ ° (aoa2 - a?) w|(q) 


(B.2) 


^L(q) = A^ + 

with the quantities a„, which are defined as 


2 Q:oa2 (w^o + Wq2) - 2a?a;2i 2 


{aoa 2 - af) 


- WG(q) 


(B.3) 


= E ((f/-)' 


H \ 


FS 


together with: 


(B.4) 


^ I] ( ^ (A^/k)' 

H \ kM 


^qn 


FS 


(B.5) 
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and the Leggett mass Ag = = 0). Thus, Eq. (IB.II) is a quadratic equation with respect to with the 

solutions 




1,2 


1 

2 


Wp(q) +w|(q) +w^(q) ± (a;|(q) + a;^(q) + ^^(q)) 


L ^ [‘^Kq)+ ^G(q)]^L (q) 

V (q) + (q) + (q)) ^ 


(B.6) 


This result can be further simplified by using a Taylor expansion of the square root. Therefore, by considering 
terms up to second order in |q| one gets: 


UJ 


2 

1,2 


= w|(q) + a;^(q) + w 2 (q) ± 


[a;|(q) +a;g;(q)] a;g(q) 
w|(q) +w^(q) +a;£(q) 


(B.7) 


In the absence of the long-range Coulomb interaction (i.e. for the case Wp(q) = 0) one finds from Eq. (IB.71) 
following result for the collective modes [see Eq. (11) in the main text]: 


ujf = ujq{c{) + 0 (—Gauge mode 
^2 = WL(q) + ^ J Leggett mode 

The Coulomb interaction leads to the renormalization of this result [see Eq. (12) in the main text]: 

ajf = ujp (q) + a;G (q) + (q) ~ (q) + ^ f T 1 Plasma mode 

/ 2/ NX (B.9) 

wf = w[2(q) + O j Leggett mode 


Thus, the mass of the Leggett mode remains unaffected by this process, but its dispersion is changed. In the 
limiting case of small q the dispersion modification is, however, negligible. 


C. ANDERSON-HIGGS MECHANISM 


In order to discuss the Anderson-Higgs mechanism for the Leggett mode in non-centrosymmetric supercon¬ 
ductors we consider the difference between the Coulomb-renormalized Leggett mode Wp defined in Eq. (IB.2|) 
and its unrenormalized counterpart wp defined in Eq. (IB.31) : 


..2.,.n ,./2.,.n «0«2 Ko + - 2Q:ja;ql - «Kl /p , N 

^Liqj iqj / 2N ^Giq] / 2 n 2 

(aoaz-af) [aoa2 - af) uj^{q) 

Eor simplicity we make following assumptions: (i) low temperature limit (T —>-0); (ii) isotropic spin-orbit 
coupling /k = 1 [corresponding to the leading order of 7k for the cubic point group 0(432)]. A generalization 
beyond these approximations is, however, straightforward. With these assumptions, equation (1C.II) simplifies to 


WL(q) - WL(q) 


2aoU;qo - 2a,uj^, ^ 

H-af) - 


2 4 2 4 

'^O^qO - 

(ag -af)a;|(q) 


(C.2) 


with Wq = i '*^F(ilqP vpfi being the Fermi velocity on the band fi = ±1. By using the definitions (IB.4|) 
and (IB.51) the equation (1C.21) can be further simplified. Thus, after straightforward calculations one obtains: 


wE(q) - wl (q) = t^G(q) 


WG(q) 


(C.3) 
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with the Fermi surface average (...) defined as (z(k'))pg = J d(f> J dOzCk') sind for a given function z(k'). Finally, 
assuming the same DoS on both bands, i.e. iV^ = Nq, and almost similar Fermi velocities, i.e. (uf+ — uf-) ^ wf 
with wf = max^ uf/^, one obtains after performing the integration: 


^L(q) -^L (q) 

u||q|2 


1 (^g+ - ^g-) (^^1+ - ^^1-) 

6 + vl_ 


(C.4) 
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